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ABSTRACT 

We discuss the role that magnetic fields in low luminosity accretion flows can play in 
creating and maintaining a multi-phase medium, and show that small magnetically- 
confined clouds or filaments of dense cold gas can dramatically reprocess the 'primary' 
radiation from tori. In particular, radio emission would be suppressed by free-free 
absorption, and an extra (weak) component would appear at optical wavelengths. 

This is expected to be a common process in various environments in the central 
regions of Active Galaxies, such as broad line regions, accretion disk coronae and jets. 
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1 INTRODUCTION 

Observational evidence, from the optical to the 7~ray band, 
suggests the presence of gas with a range of physical 
properties in the inner regions of Active Galaxies (AGN). 
Hot/relativistic material co-exists with 'cold' quasi-thermal 
gas, on a large range of scales in these complex environ- 
ments: tenuous hot plasma forms a corona above or around 
relatively cold and dense plasma, in the form of an accretion 
disk or re-processing clouds; clumps of photoionized gas, re- 
sponsible for the observed narrow and broad emission lines, 
are embedded in a tenuous hotter medium. 

The formation, confinement and survival of a multi- 
phase medium has been considered with much detail in the 
case of thermal re-processing by the material which gives 
rise to the broad lines. Difficulties associated with the main- 
tenance of photoionized matter in two phases which are in 
gas pressure equilibrium, first led to the suggestion that a 
more effective way of sustaining temperature and density 
gradients may involve magnetic fields in equipartition with 
other forms of energy (Rees 1987). 

Magnetic fields are believed to be a common ingredi- 
ent in the AGN environment and so can provide effective 
confinement and insulation of cold gas embedded in a much 
hotter ambient medium. In particular, more detailed work, 
focused on smaller scales, has shown that this is plausibly 
the case in the magneto-sphere around the central engine 
(typically on scales < 50i? s ) as well as in relativistic jets, 
thanks to the combined action of the magnetic and radia- 
tion fields (Celotti, Fabian & Rees 1992; Kuncic, Blackman 
& Rees 1996; Kuncic, Celotti & Rees 1997; Celotti et al. 
1998; Celotti & Rees 1998). 

Recently, much attention has focussed on accretion 



flows where electron-ion coupling could be ineffective be- 
cause the density is low and the kinetic temperature very 
high. Very hot ions then provide the dominant pressure, but 
radiative cooling is inefficient because the energy cannot be 
tapped by the (cooler) electrons (Rees et al. 1982). Accre- 
tion flows of this kind, in which thermal energy is advected 
inwards rather than being radiated, may be relevant to the 
puzzling absence of substantial radiative dissipation in large 
elliptical galaxies (Fabian & Rees 1995) . This has raised in- 
terest in the phenomenon and led to more detailed theo- 
retical work and the comparison of the expected spectral 
signatures with data. Approximate analytical solutions as 
well as self-consistent numerical results for the structure and 
properties of the flow (re-named ADAF) have been derived 
(e.g. Narayan & Yi 1994, 1995; Chen et al. 1995; Mahadevan 
1997; Chen, Abramowicz & Lasota 1997). Although the faint 
ADAF emission cannot be clearly disentangled from other 
spectral components over a significant energy range, the ra- 
diative predictions seem to be, so far, globally consistent 
with data (e.g. Lasota et al. 1996; Reynolds et al. 1996; Di 
Matteo & Fabian 1997; Esin, McClintock & Narayan 1997; 
Narayan et al. 1998). However, thanks to very recent ob- 
servations of a few elliptical galaxies in the radio to mm 
frequency range, Herrnstein et al. (1998) and Di Matteo et 
al. (1998) have pointed out that (at least) 'standard' ADAF 
models might be inadequate to account for the (absence of) 
emission in large elliptical galaxies. 

Here we suggest the possibility that within the flow 
some of the gas can be in the form of cold and dense clumps 
or filaments, confined by the magnetic field expected to 
thread the accreting matter and examine the effects that 
this might produce on the observed spectral signatures of 
ADAFs. 
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2 THE MODEL AND ITS INGREDIENTS 

Let us consider the properties of an ADAF, according to 
the accretion solutions derived by Narayan & Yi (1995) and 
Mahadevan (1997). The self-similarity postulated by these 
authors must break down near the inner boundary, but may 
nevertheless be an adequate approximation over a substan- 
tial range in r. 

The expected field, assumed in quasi-equipartition with 
the total gas pressure, is of the order of 

B = 3.7 x lO'm-^m^V- 5 / 4 G (1) 

where mg = M/10 s Mq is the mass of the black hole, 
m_ a = M/1CT 3 Me the accretion rate in Eddington units 
and rs = R/5Rs the distance in units of Schwarzschild radii. 
Hereafter, the viscosity parameter and the gas to total pres- 
sure have been taken as a ~ 0.3 and j3 — 0.5, respectively, 
though it is straightforward to scale our results to other val- 
ues (the constants ci and C3 of the order of unity as defined 
in Narayan & Yi 1995, are taken ci = 0.5 and C3 = 0.3). 

2.1 Dense and cold gas 

The magnetic field can trap some plasma and, as this is 
likely to have cooled to the Compton temperature, a density 
contrast as high as ~ T v i r /Tc m P t with respect to the exter- 
nal/hot phase medium can be attained. Two body processes 
start to become very efficient in cooling the dense trapped 
gas, which can thus cool still further, typically reaching tem- 
peratures of the order of T ~ 10 4 24 K (see Section 3.2). 

The gas pressure is limited by the balance with the total 
magnetic field stresses and therefore a maximum density for 
the cold phase, achieved for a diamagnetic clump of gas, can 
be as high as 

n ~ n hot T vir /T ~2x 10 16 m^ 1 m- 3 r~ 3/2 T^ 1 cm" 3 , (2) 
where n bot represents the density in the hot phase. 

2.2 Geometrical structure 

The spatial distribution of the gas is clearly determined by 
the field structure and the net external force acting on it. 
This would limit the scale-height of the gas to a typical 
value such that the tension of the distorted field (in the 
perpendicular direction) balances the net effect of gravity, 
radiation pressure and acceleration. Typically, in presence of 
loops with a toroidal field component, the action of gravity 
leads to dimensions h g ~ kT/m p g. However, other forces 
can dominate in these environments, such as the gas inertia 
to dynamical acceleration and the radiation force. Despite 
the significant opacity (mostly free-free) of the extremely 
dense and cold gas, because the low ADAF luminosity we 
do not expect radiative acceleration to play a major role in 
this environment (see Section 3.2), in contrast with other 
situations in the central regions of AGN. 

It should be stressed that, from the point of view of 
the large-scale dynamics, these structures of cold and dense 
gas can be treated as a uniform loading to the global flow, 
as long as the individual clouds/filaments are small enough 
that, dragged by the field, they are constrained to follow the 
same dynamics as the hot material. This requires them to be 
smaller than the scale-height h e g imposed by the effective 



local acceleration a, which may be dominated by the effect 
of magnetic stresses or internal random motions (see also 
Section 3.1). As these are difficult to quantify, we have a 
fiducial upper limit to the typical scale-height of the order 
of ft g 

h e g ~ 10 6 -m<T 1 Tat\ cm. (3) 
a 

3 THE SURVIVAL OF THE COLD PHASE 

Let us consider the survival of clumps of cold matter thus 
confined and advected with the flow. 

3.1 Diffusive processes in clouds and filaments 

At the high densities inferred, radiative (see Section 3.2) and 
collisional timescales are much shorter that the dynamical 
one, estimated as the free fall time, iff ~ 10 4 mg ll>2 rfJ 2 s, 
allowing the gas to reach thermal and radiative equilibrium. 
Also, pressure equilibrium with the surrounding medium can 
be easily maintained in the global dynamical timescale, for 
typical dimension of the cold phase structures smaller than 

< inlO -l/2rpl/2 3/2 

£ 10 m g ' T 4 ' r s ' cm. 

Although h e « is much smaller that the extension of the 
hot gas, it greatly exceeds, for the inferred field strength, 
the characteristic gyro-radii of particles in both phases. 

Resistivity is low enough to maintain the coupling be- 
tween field and gas on radiative and dynamical timescales. 
However, collisions of the hot particles with the dense 
clumped gas, can lead to effective thermal conduction. As 
a detailed radiative computation has shown (see Section 
3.2), the cold gas is only marginally ionized and, under 
this circumstance, the dominant diffusion is due - at least 
in the inner ADAF regions - to interactions with neutral 
atoms, with a typical mean free path A ~ 0.1nj~ 6 cm. In a 
timescale i the cloud depth affected by thermal conduction, 
according to the diffusion approximation, is of the order of 
d ~ 4 x 10 4 n^ 6 1/ ' 2 i 1 ' /2 cm. Furthermore, the collision fre- 
quency is so high with respect to the gyration frequency 
(typically v C oii ~ 10 fgyr) that even in the direction per- 
pendicular to the field line, conduction is effective (the col- 
lisional frequency is proportionately lower on larger scales). 
Even when clouds are thinner than the diffusion length, ra- 
diative cooling may allow them to survive. 

In a disc corona or a relativistic jet the field can be ei- 
ther anchored in the disk (corona) or have an ordered com- 
ponent (jet). A possibly severe threat to filament survival 
- as already mentioned - is posed in an environment where 
random motions of field loops and large-amplitude Alfven 
waves would cause violent local accelerations, which corre- 
spondingly reduce the thickness of filaments that could be 
carried with the flow. The acceleration would be largest on 
the smallest scales. Let us consider, simply as an illustra- 
tion, a Kolmogorov spectrum where the velocities associated 
with a scale I are proportional to l 1 ^ 3 : the associated accel- 
erations would then depend on scale as l^ 1 ^ 3 . This means 
that the filaments of a given scale can less easily follow the 
smaller-scale eddies. On the other hand, stretching of field 
lines would decrease the scale of the filaments. 

We draw attention to this issue, emphasizing the un- 
certainty in the fate of dense filaments that are advected in- 
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wards with the flow towards regions of progressively stronger 
field nearer the hole. One can envisage an element of gas, 
trapped within a magnetic tube, starting off in the broad 
line region and moving into the denser and more extreme 
conditions envisaged in the ADAF. But can the cool-phase 
gas survive, or will it be shredded and homogenized? It is 
unclear whether the disordered motions associated with the 
magnetic viscosity would shred any filaments into structures 
so thin that they cannot survive the effects of conductivity 
and diffusion. 

This is an important open question for matter which is 
coupled to the field lines. However, material decoupled from 
the field (e.g. ejecta from stars, behaving like diamagnetic 
clouds) could survive more easily since it could be squeezed 
by the surrounding field but would not necessarily follow 
the violent random motion associated with the 'a-viscosity'. 
Therefore, irrespective of the answer to the question posed 
above, small clouds or filaments of magnetically-confined 
cool gas could still have important consequences for the 
ADAF spectrum, as shown in the next Section. 

3.2 Radiative signatures 

Because of the tiny transverse dimensions, the cold matter 
occupies a small fraction of the total volume, even though 
a significant amount of the total mass inflow might re- 
side in this component if its scale-height is comparable 
to h g (~ (T/T v i r )R). Nevertheless, because of the extreme 
high density, this matter can be very effective in absorb- 
ing (mainly through free-free) and re-processing radiation 
impinging on it. 

Furthermore, the 3D structure of the field can allow the 
matter to accumulate on sheet-like surfaces or filaments in 
the direction perpendicular to the net force, distributing it 
with high covering factors. Even though the volume filling 
factor is small, the surface covering factor can easily be of 
order unity. 

Because of the comparatively low radiation field in 
ADAFs, the gas might cool to temperatures of the order of 
10 4 K at ~ 57?s, allowing a density contrast between the two 
phases as high as 10 8 . The low temperature and high den- 
sity cause the trapped gas to be an efficient absorber, even 
more so as a significant fraction of the ADAF radiation is 
emitted as thermal cyclo-synchrotron at low frequencies. 

Therefore, if dense filaments existed within an ADAF, 
with (at least at some radii) a substantial covering factor, 
they could modify the standard ADAF spectrum. In partic- 
ular, the radio and mm bump could be absorbed by free-free 
(re-appearing in the optical or UV bands). [Furthermore, for 
sufficiently ionized gas, the cold phase plasma frequency also 
exceeds that of the cyclo-synchrotron peak, raising the pos- 
sibility of collective absorption or reflection effects.] The ac- 
tual result on the radio spectrum would depend on whether 
the clouds were just in the inner parts (where they would ab- 
sorb only the higher-frequency sub-mm radiation) or mainly 
in the outer part (in which case they could absorb/reflect 
everything), where the gas would be largely molecular and 
grain could survive. We stress here that this is an interest- 
ing possible complication in confronting ADAF models with 
observations. 

Under these circumstances, the absorbed radiation can 
be energetically important and the spectral effect rele- 
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Figure 1. Reprocessed spectra emitted by cold gas structures lo- 
cated in an ADAF environment at about 5-Rg and with a covering 
factor ~ 1. The dashed and continuous lines indicate the incident 
and transmitted spectra, respectively. 



vant. We computed the effective temperature and conditions 
reached by the cold gas in thermal and radiative equilib- 
rium, as well as the reprocessed spectrum from a detailed 
radiative transfer calculation, by using the code CLOUDY 
(it should be however noted that the code results might be 
not totally reliable in these extreme conditions). The local 
incident spectrum has been assumed to be that of an ADAF 
at the considered distance (e.g. Mahadevan 1997), including 
the synchrotron, Compton and bremsstrahlung components. 

Even in the inner part of the ADAF (~ 5Rs) the equilib- 
rium temperature, as already mentioned, only reaches ~ 10 4 
K and the gas is largely in neutral form. 

In Fig. 1 we show the result of the computation when 
a covering factor of ~ 1 is assumed. Substantially all of the 
cloud heating is due to free-free absorption of the cyclo- 
synchrotron component, which is balanced by collisional 
line and free-bound cooling, the longest radiative timescale 
being ~ 50 s. The spectral features in the re-processed 
spectrum are expected to be broadened and shifted by the 
high velocity motion of the gas. This includes any (weak) 
iron emission line, which therefore would be probably unde- 
tectable. The column density in neutral hydrogen can be less 
than 10 22 cm -2 , thus affecting only the soft X-ray emission 
from the ADAF. 

Although the ionization state and the spectral results 
clearly depends on the detailed incident spectrum (e.g. the 
amount of Comptonized synchrotron photons), the result 
shown here is representative of the global effect of the re- 
processing by the cold material. 
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4 DISCUSSION AND CONCLUSIONS 

We have stressed the key role that magnetic fields in quasi- 
cquipartition in the accretion flow of AGN can have in cre- 
ating and maintaining a multi-phase medium, by confining 
cold, dense gas with small filling factor. This is expected to 
be a common process in the central regions of Active Galax- 
ies (broad line region, accretion disk coronae, jets,...). 

As mentioned in the introduction, high energy observa- 
tions of AGN strongly suggest the presence of gas in different 
physical conditions in the very center of AGN. In particu- 
lar, spectral evidence has accumulated for a component of 
optically thick, cold gas, as well as a tenuous hot medium, 
to account for the ultraviolet to 7-ray spectrum of radio- 
quiet AGN and galactic black hole candidates. Interestingly, 
Krolik (1998) has recently proposed that a consistent and 
stable solution for accretion disc might involve the formation 
of clumps of dense gas -albeit warmer and optically thicker 
than those considered here- embedded in a hot medium and 
threaded by magnetic field lines. 

We have shown that thermal gas trapped by the field in 
an ADAF, can survive in the hot environment long enough 
to significantly reprocess the 'primary' radiation and pro- 
duce observable spectral signatures. Its predominant effect 
would be to largely suppress the cyclo-synchrotron compo- 
nent produced by the hot gas, in the radio-mm band, which 
has been considered, so far, an observable indication of the 
presence of radiative inefficient accretion phase in our Galac- 
tic Centre and in other low-luminosity galactic nuclei (and 
Galactic Black Hole candidates). The resulting spectrum 
clearly depends on the radial distribution of the amount and 
covering factor of the gas. There are here two countervailing 
effects that are hard to quantify. The pressures (and hence 
the maximum densities of cool gas) are higher at small radii, 
so a given amount of gas would cause more absorption near 
the centre. On the other hand, if the cool gas consists of fil- 
aments entrained in the flow, and advected inwards with it 
(rather than being free of magnetic flux and merely squeezed 
and confined by the hot-phase medium) , the violent random 
motions associated with the magnetic viscosity may shred 
and disperse it before it reaches the centre: in this case the 
dominant effect would be due to gas at 100-1000 Rs- 

Interestingly, the presence of cold clumps might account 
for the tight upper limits recently found on the radio, cm 
and sub-mm emission of large elliptical galaxies likely to host 
a massive black hole accreting ISM in the radiatively inef- 
ficient ADAF regime. In particular the observational con- 
straints on the 22 GHz emission from the region within 
the molecular maser disk in NGC 4258 (Herrnstein et al. 
1998) have been strengthened by the very recent results 
of VLA and SCUBA observations of three elliptical galax- 
ies by Di Matteo et al. (1998). The upper limits found by 
these authors are indeed well below the expected emission 
from 'standard' ADAF models: cold material in the condi- 
tions considered here could be responsible for the absorption 
of a significant amount of the (unobserved) ADAF cyclo- 
synchrotron radiation. 

The strong lines and edges, but especially the thermal 
continuum shape in the optical-UV band, expected from the 
reprocessing by cold material, might be soon observable, as 
demonstrated by the upper limits already derived from the 
scattered optical light in NGC 4258 (Wilkes et al. 1995). 
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